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Abstract 

The two-point function and the vacuum expectation value (VEV) of the current density are 
investigated for a massive charged scalar field with arbitrary curvature coupling in the geometry of 
a brane on the background of AdS spacetime with partial toroidal compactification. The presence of 
a gauge field flux, enclosed by compact dimensions, is assumed. On the brane the field obeys Robin 
boundary condition and along compact dimensions periodicity conditions with general phases are 
imposed. There is a range in the space of the values for the coefficient in the boundary condition 
where the Poincare vacuum is unstable. This range depends on the location of the brane and is 
different for the regions between the brane and AdS boundary and between the brane and the 
horizon. In models with compact dimensions the stability condition is less restrictive than that for 
the AdS bulk with trivial topology. The vacuum charge density and the components of the current 
along non-compact dimensions vanish. The VEV of the current density along compact dimensions 
is a periodic function of the gauge field flux with the period equal to the flux quantum. It is 
decomposed into the boundary-free and brane-induced contributions. The asymptotic behavior of 
the latter is investigated near the brane, near the AdS boundary and near the horizon. It is shown 
that, in contrast to the VEVs of the field squared and energy-momentum tensor, the current density 
is finite on the brane and vanishes for the special case of Dirichlet boundary condition. Both the 
boundary-free and brane-induced contributions vanish on the AdS boundary. The brane-induced 
contribution vanishes on the horizon and for points near the horizon the current is dominated by 
the boundary-free part. In the near-horizon limit, the latter is connected to the corresponding 
quantity for a massless field in the Minkowski bulk by a simple conformal relation. Depending on 
the value of the Robin coefficient, the presence of the brane can either increase or decrease the 
vacuum currents. Applications are given for a higher-dimensional version of the Randall-Sundrum 
1-brane model. 


PACS numbers: 04.62.+V, 04.50.-h, ll.10.Kk, ll.25.-w 


1 Introduction 

Anti-de Sitter (AdS) spacetime is one of the simplest and most interesting spacetimes allowed by 
general relativity. It is the unique maximally symmetric solution of the vacuum Einstein equations 
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with a negative cosmological constant (for geometrical properties of AdS space and its uses see, 
e-g-, IB)- Quantum field theory in AdS background has long been an active field of research for a 
variety of reasons. First of all, AdS spacetime has a high degree of symmetry and, because of this, 
numerous physical problems are exactly solvable in this geometry. The maximal symmetry of AdS 
simplifies many aspects of the study of quantum fields and the investigations of the corresponding field- 
theoretical effects may help to develop the research tools and insights to deal with more complicated 
geometries. Much of early interest to quantum field theory on AdS bulk was motivated by principal 
questions of the quantization of fields on curved backgrounds. The lack of global hyperbolicity and 
the presence of both regular and irregular modes give rise to a number of new features which have 
no analogues in quantum held theory on the Minkowski bulk. Being a constant negative curvature 
manifold, AdS space provides a convenient infrared regulator in interacting quantum held theories [2]. 
Its natural length scale can be used to regularize infrared divergences without reducing the symmetries. 
The importance of this theoretical research was increased by the natural appearance of AdS spacetime 
as a ground state in supergravity and Kaluza-Klein theories and also as the near horizon geometry of 
the extremal black holes and domain walls. 

A further increase of interest is related to the crucial role of the AdS geometry in two exciting 
developments of the past decade such as the AdS/CFT correspondence and the braneworld scenario 
with large extra dimensions. The AdS/CFT correspondence [3] (see also 0) represents a realization of 
the holographic principle and relates string theories or supergravity in the AdS bulk with a conformal 
field theory living on its boundary. It has many interesting consequences and provides a powerful 
tool for the investigation of gauge held theories. Among the recent developments of the AdS/CFT 
correspondence is the application to strong-coupling problems in condensed matter physics (familiar 
examples include holographic superconductors, quantum phase transitions and topological insulators). 
Moreover, the correspondence between the theories on AdS and Minkowski bulks may be used to derive 
new results in mathematical physics, in particular, in the theory of special functions (see, for instance, 
[5] and references therein). The braneworld scenario (for reviews see |6]) offers a new perspective on 
the hierarchy problem between the gravitational and electroweak mass scales. The main idea to resolve 
the large hierarchy is that the small coupling of four-dimensional gravity is generated by the large 
physical volume of extra dimensions. Braneworlds naturally appear in string/M-theory context and 
present intriguing possibilities to solve or to address from a different point of view various problems 
in particle physics and cosmology. 

An inherent feature of all these models is that the boundary conditions on the fields should be 
specified in order to completely determine the dynamics. First of all the AdS spacetime is not globally 
hyperbolic and has a time-like future null infinity. As a consequence of this, the information may be 
lost to, or gained from, spatial infinity in finite coordinate time. In order to define a consistent quantum 
field theory, appropriate boundary conditions must be imposed BIB!. The general class of allowed 
boundary conditions on the AdS boundary has been discussed in [9], based on the analysis of Hi- 
Different boundary conditions lead to different theories, in particular to different conformal field theory 
duals. In braneworld models on AdS bulk, the presence of branes gives rise to additional boundary 
conditions on the operator of a quantum held. These conditions depend on the specific geometry of 
the brane and have been discussed in mmm for models with Z 2 symmetry. Another type of 
boundary conditions is induced by the presence of compact spatial dimensions. The extra compact 
dimensions are an inherent feature of braneworld models arising from string and M-theories. In 
these models one needs also to specify the periodicity conditions along compact dimensions. Different 
conditions correspond to topologically inequivalent held configurations (14j . The nontrivial periodicity 
conditions lead to a number of interesting quantum held-theoretical effects, which include instabilities 
in interacting held theories, topological mass generation, and symmetry breaking. 

All these types of boundary conditions, imposed on the held operator, modify the spectrum of the 
zero-point fluctuations of a quantum held. As a consequence, the vacuum expectation values (VEVs) 
of physical observables are shifted by an amount depending on the geometry of the boundary and on 
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the type of the boundary condition. This general phenomena is known as the Casimir effect. It has 
important implications on all scales, from mesoscopic condensed matter physics to cosmology. Since 
the original research by Casimir (15j . for the electromagnetic field in the geometry of two parallel 
conducting plates, many theoretical and experimental works have been done in this direction (for 
reviews see [16]). In braneworlds, the boundary conditions imposed on the bulk fields will give Casimir- 
type contributions to the vacuum energy and to the vacuum forces acting on the branes. The latter 
provide a natural mechanism for stabilizing the intei'brane distance (radion) in Randall-Sundrum-type 
models. The Casimir energy gives a contribution to both the brane and the bulk cosmological constants 
and should be taken into account in the self-consistent formulation of the scenario. Motivated by these 
issues, the investigations of the Casimir energy and related forces on AdS bulk have attracted a great 
deal of attention (see, for instance, the references in B2D- The Casimir effect in higher-dimensional 
generalizations of the AdS spacetime with compact internal spaces has been discussed in [18) . 

An important physical characteristic of the vacuum state for charged fields is the expectation value 
of the current density. It carries information about the geometry and topology of the background space 
and is responsible for the backreaction of the quantum field, as a source in semiclassical Maxwell’s 
equations. In the present paper we investigate the VEV of the current density for a charged scalar 
field in background of locally AdS spacetime with an arbitrary number of toroidally compactified 
spatial dimensions, in the presence of a brane parallel to the AdS boundary. The corresponding 
problem in the absence of the brane has been considered recently in [19] and here we shall be mainly 
concerned with the brane-induced effects. Both the zero and finite temperature expectation values 
of the current density for charged scalar and fermionic fields in background of flat spacetime with 
toroidal dimensions were investigated in muzE]. Applications were given to the electronic subsystem 
of cylindrical and toroidal carbon nanotubes described in terms of a (2 + l)-dimensional effective field 
theory. The vacuum current densities for charged scalar and Dirac spinor fields in de Sitter spacetime 
with toroidally compact spatial dimensions are considered in [22], The influence of boundaries on 
the vacuum currents in topologically nontrivial spaces are studied in [23, 23 for scalar and fermionic 
fields. The effects of nontrivial topology induced by the compactification of a cosmic string along its 
axis have been discussed in [25]. 

The organization of the paper is as follows. In the next section we specify the bulk and boundary 
geometries under consideration and evaluate the Hadamard function for a charged massive scalar 
field in both regions on the right and on the left to the brane (referred to as R- and L-regions, 
respectively). The brane-induced contributions are manifestly extracted and they are presented in 
the form well suited for the investigation of the VEVs of local physical observables bilinear in the 
field. As such an observable, in section [3[ we consider the current density. The corresponding VEVs 
are decomposed into boundary-free and brane-induced contributions for both R- and L-regions. The 
asymptotic behavior of the brane-induced contributions is considered near the brane, near the AdS 
boundary and near the horizon. Limiting expressions are derived for small and large proper lengths 
of compact dimensions. The main results are summarized in section [4] In Appendix we provide 
alternative representations for the Hadamard functions in the R- and L-regions. The expressions for 
the current densities obtained from these representations are used for the investigation of the near 
brane asymptotic. 

2 Geometry of the problem and two-point functions 

2.1 Set-up 

We consider a charged quantum scalar field <p(x) with the mass rn and with the curvature coupling 
parameter £. In the presence of an external classical gauge field , the corresponding field equation 
reads 

+ m 2 + £R) ip(x) = 0, (2.1) 
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where = V^+ieA^, e is the charge of the field quanta, is the operator of the covariant derivative 
associated with the metric tensor g^, and R is the Ricci scalar for the background spacetime. The 
background geometry in the present paper is given by the interval 

ds 2 = g^dx^dx" = e~ 2y 1 ' a rjikdx 1 dx k — dy 2 , (2.2) 

where 77 ^ = diag(l, —1,... , —1) is the metric tensor for D-dimensional Minkowski spacetime, i,k = 
0,..., D — 1, and y, v run from 0 to D. The local geometry described by (12.21) coincides with that 
for (D + l)-dimensional AdS spacetime of the radius a, expressed in Poincare coordinates. The corre¬ 
sponding metric tensor is a solution of the Einstein equations with a negative cosmological constant 
A = —D(D — l)a~ 2 /2 and for the Ricci scalar one has R = —D(D + 1 )/a 2 . 

The spatial topology considered here will be different from that for AdS. Namely, we assume that 
the subspace with the coordinates x l , l = p + 1,..., D — 1, is compactified to a (/-dimensional torus 
T q with q = D — p— 1 (for a recent review of quantum field-theoretical effects in toroidal topology see 
[26]). We shall denote by Li the length of the Zth compact dimension, 0 ^ x l ^ L[. The ranges of the 
remaining coordinates are —00 < x l < + 00 , l = 1, 2 ,... , p, and —00 < y < + 00 . Consequently, for the 
subspace perpendicular to the y-axis we take the topology RP x T q . Note that Li is the coordinate 
length of the compact dimension. For a fixed y, the proper length of the 1th compact dimension is 
given by = e~ y ^ a Li and it decreases with increasing y. The coordinate transformation 

z = ae y / a , 0 ^ z < 00 , (2-3) 

brings the interval () 2 . 2 |) into manifestly conformally flat form with the conformal factor ( a/z ) 2 : 

ds 2 = (a/z) 2 (r]ikdx l dx k — dz 2 ). (2.4) 

In terms of the new coordinate z, the AdS boundary and horizon are presented by the hypersurfaces 
z = 0 and z = 00 , respectively. For the proper length, measured by an observer with a fixed coordinate 
z, one gets Lorn = aLi/z. 

We consider a field theory in a non-simply connected spacetime and for the complete formulation 
of the problem the periodicity conditions along compact dimensions should be specified. Here we 
impose the conditions 


<f(t, x 1 ,..., x l + L u ..., y) = e iai tp{t , x 1 ,..., x 1 ,..., y), (2.5) 

for l = p + 1,D — 1 with constant phases og. In the literature, the most frequently considered 
special cases correspond to ag = 0 (untwisted scalar) and og = 7 r (twisted scalar). As we will see below, 
the nontrivial phases in the periodicity conditions give rise to the vacuum currents along compact 
dimensions. In the discussion below we shall assume that the gauge field A^ is constant. Though the 
corresponding held strength vanishes, because of the nontrivial topology of the background spacetime, 
the VEVs of physical observables will be influenced by this sort of held configuration. This is an 
Aharonov-Bohm like effect of a constant gauge held. Let us consider two sets of the helds (ip(x),A^) 
and (ip'(x), A' ) connected by the gauge transformation <p(x) = e~ tex ^ x >ip'{x), A M = A' fl + c^x(x). For 
a constant gauge held, taking the function x( x ) = we see that in the new gauge the vector 

potential vanishes, A' fI = 0. But, after the gauge transformation the vector potential of the former 
gauge appears in the periodicity conditions for the new held operator: 

ip'it, s 1 ,..., x l + L h ..., y) = e*°V(f, x 1 ,..., x l ,..., y), ( 2 . 6 ) 

with l = p + 1,... ,D — 1, and with the new phases 


og = og + eAiLi. 


(2.7) 
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In particular, nontrivial phases are generated for untwisted and twisted scalars. The phase shift in 
(12.61) is related to the magnetic flux T; enclosed by the Ith compact dimension: eA[Li = 
where d>o = 2-7T/e is the flux quantum. Note that the gauge field fluxes play an important role in 
recent developments of string theory compactifications (for a review see [27]). 

In the problem under consideration, the second type of boundary condition imposed on the field 
operator is induced by a brane parallel to the AdS boundary and located at y = r/o- The corresponding 
value for the conformal coordinate z we shall denote by zo = ae yo ^ a . On the brane we assume a gauge 
invariant boundary condition of the Robin type: 

(1 + = 0, y = y 0 , (2.8) 

where /3 is a constant and n y is the inward pointing normal to the brane. For the latter one has 
n y = in the region y > yo and n y = —5^, in the region y < yo- Note that, in general, the value of the 
coefficient j3 for these regions could be different. The Robin boundary condition is a generalization of 
Dirichlet and Neumann conditions and naturally appears in a number of physical problems, including 
those in braneworld scenario (see below). The spatial geometry under consideration for D = 2, 
embedded into the 3-dinrensional Euclidean space, is displayed in figure [1] The compact dimension 
is presented by the circles and the thick circle corresponds to the location of the brane. We have 
also depicted the gauge field flux tube enclosed by the compact dimension. The proper length of the 
compact dimension decreases with increasing y. 



Figure 1: The spatial section of the geometry for D = 2 embedded into a 3-dimensional Euclidean 
space. The thick circle corresponds to the brane. 

The physical quantity we are interested in is the VEV of the current density, (0|j /i (x)|0) = (j M (x)), 
for the field ^p(x), where |0) stands for the vacuum state. In quantum field theory on curved back¬ 
grounds the choice of the vacuum state is among the basic points. In what follows it will be assumed 
that the field is prepared in the Poincare vacuum state. The latter is realized by the mode functions 
of the field which are obtained by solving the field equation in Poincare coordinates (for the discussion 
of the relation between the Poincare and global vacua see, for instance, [28]). The operator of the 
current density for a charged scalar field is defined by the expression 

3n{x) = ie[ip + (x)D^(p(x) - {D(x))ip(x)}. (2.9) 

Its VEV is among the most important characteristics of the vacuum state. The procedure we shall 
use here for the evaluation of the expectation value is based on the formula 

(jfAx)) = lim (d M - d' + 2 ieA^)G{x, x'), (2.10) 

Z x'^tx 

where 

G(x,x' ) = (0|(/?(x)<£ + (a/) + v9 + (x , )<y9(x)|0) 

= EE^’wV’V)- (2.n) 

a s=± 
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is the Hadamard function for the vacuum state under consideration. In (12. lip , the summation goes 
over a complete orthonormal set of positive- and negative-energy mode functions {x) specified 
by a collective index a involving the corresponding quantum numbers. The mode functions obey the 
quasiperiodicity conditions (12.51) and the condition (12.81) on the brane. 

So, as the first step, we shall evaluate the Hadamard function. Though the background AdS 
spacetime is homogeneous, the brane at y = yo has nonzero extrinsic curvature tensor and its sides are 
not equivalent. In particular, the VEVs differ in the regions on the right and on the left of the brane. 
The consideration requires different procedures for these regions and we discuss them separately. In 
what follows the regions y > yo and y < yo will be referred to as R-region (right region) and L-region 
(left region), respectively. 

2.2 Hadamard function in the R-region 

We shall work in the gauge with the fields {tp'(x), A'^ = 0), omitting the prime. In accordance with 
the problem symmetry, the mode functions can be factorized as 

^±\x) = z D/2 Z u {Xz)e ikrXrTiut , ( 2 . 12 ) 

where the summation over r in the exponent goes for r = 1,... ,D — 1, Z u (x) is a cylinder function of 


the order 

z, = y/02/4 - D(D + 1)£ + m 2 a 2 , 

(2.13) 

and 

D -1 



uj = VA 2 + k 2 , k 2 = kf. 

i=i 

(2.14) 


For imaginary values of the order v the vacuum state becomes unstable [8j and in what follows we 
shall assume the values of the parameters for which v ^ 0. In the cases of conformally and minimally 
coupled massless fields one has v = 1/2 and v = D/2, respectively. In the former case, we have 
the standard conformal relation with the modes in the problem on Minkowski bulk with toroidal 
dimensions. In (|2.12|) . for the components of the momentum one has —oo < hi < +oo, l = 1,... ,p, 
and 

h = (2irni + &i)/Li, l = p + 1,..., D - 1, (2.15) 

with rii = 0, ±1,±2,.... The eigenvalues (12.151) for the components along compact dimensions are 
directly obtained from the quasiperiodicity conditions (|2.6I) . 

In the R-region, the function Z v {Xz) is a linear combination of the Bessel and Neumann functions, 
J u (Xz) and Y u (Xz). First we consider the case when for all the modes A is real. The changes in the 
evaluation procedure in the case when the modes with purely imaginary values of A are allowed (bound 
states) will be discussed below. The relative coefficient in the combination of the Bessel and Neumann 
functions is determined from the boundary condition (12.81) on the brane (with D^ = d ] /l ) and one gets 

Z v {Xz) = C a g v (XzQ, Xz), (2.16) 

where, for the further convenience, we have introduced the function 

gu(u,v) = J u (v)Y v (u) - J v (u)Y v (v). (2.17) 

Here, for a given function F{x ), we use the notation 

F{x) = AqF{x) + BqxF'(x), (2.18) 

with the coefficients 

A 0 = 1 + y 5 y P/a , B 0 = 5 y f3/a. (2.19) 
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In (|2.19l) and in what follows, 5 y = 1 in the R-region and S y = — 1 in the L-region. 

Now, the set of quantum numbers a is specified by a = (A, k p , n,j), where k p = {k \,..., k p ) is the 
momentum in non-compact space and rig = (n p + 1 ,... ,riD-i)- The coefficient C a is determined from 
the normalization condition 


d D xg 00 y/\g\<p { j\x)v i ° ] *(x) = -^<5(A - A')<5(k p - k(Jh r 


L q 7 


( 2 . 20 ) 


where the y-integration goes over [yo, oo). By using the mode functions (12.121) with the radial function 
from (|2.16l) . we find 

l<?a| 2 = 2^^ [ J Z( Xz o) + Y?(\zo)]-\ (2.21) 

with V q = L p+ 1 • • • Ld -i being the volume of the compact subspace. 

Having determined the complete set of normalized modes, from the mode-sum in (12.111) . for the 
Hadamard function one obtains the representation 


G(x,x') = 


(zz') 


l\D/2 


£ 


dk p e' 


ik r Ax r 


dXX 


(2vr ) p a D -v q — J Jo 

n q 

_ cos(AtVA 2 + k 2 ) g u (Xz 0 , Xz)g u (Xz 0 , Xz’) 
YW+W J 2 (Xz 0 ) + Y 2 (Xz 0 ) 


( 2 . 22 ) 


In order to extract explicitly the brane-induced contribution, we subtract from (12.221) the correspond¬ 
ing Hadamard function in the boundary-free AdS background which will be denoted by Gq(x,x'). 
The latter is obtained from (12.221) replacing the last fraction in the right-hand side by the product 
J u (Xz)J u (Xz'). By using the identity 


- 1 £ A ))*’ 1 

f2/w„)iy2(wA 2 ^ 2 h1 3) {Xz 0 ) 


tf{Xzo) + i; 2 (Az 0 ) 

/1 2 ^ 

with Hy ’ ; ( x ) being the Hankel functions, the following decomposition is obtained 

/ roo 

dk p e ikrAxr j dX 


a 1 - 0 (zz') D/2 

G(x,x') = Go{x,x') — 


2 (27 T) P V q 


','CPS(A WX + k?) y U\z a ) HC ,) (Az)gfa)(A/) 

Xzo) 


(2.23) 


(2.24) 


Now, assuming that z > zo, we rotate the integration contour of the last integral by the angle n/2 
for j = 1 and by — tt/2 for j = 2. Introducing the modified Bessel functions I v (x) and K u (x), the 
following final expression is obtained: 


4 ( p 

G(x,x') = GoM - j dk 

X / 


Ak r Ax r 


°° cosh(Af y/u 2 - k 2 ) I u {uzo) T ,. < , ^ 

duu - , -- 7 I\ v [uz)K u [uz ), 


Vu 2 - k 2 Ky(uzo)' 


(2.25) 


where the notation with overbar is defined by (12.181) . (12.191) with S y = 1. 

In deriving (|2.25p we have assumed that for all the modes of the field A is real. In addition to 
these modes, bound states can be present. For them A is purely imaginary, A = ig, g > 0, and the 
mode functions have the form 


<p^{x) = C { ^z D l 2 K v {j ] zY KxT ^ iu( ^\ 


(2.26) 
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where oj{rj) = a / k 2 — r / 2 . Let us denote by k^j the lowest value of the momentum in the compact 
subspace. Assuming that |cb| ^ ir, one has 

k (°J 2 = E ( 2 - 27 ) 

i=p +1 

If r] > k^j, then there are modes for which the energy is purely imaginary and the vacuum state 

is unstable. In order to have a stable vacuum, in what follows we assume that r/ < k^j. From the 
boundary condition (12.81) it follows that for bound states the possible values of rj are roots of the 
equation 

K u (vzo) = 0. (2.28) 

By using the recurrence relation for the Macdonald function, this equation can be rewritten as 


{v - D/2 — a/)3) K u (u) + uK u _i(u) = 0, 


(2.29) 


with u = tjzq. From here it follows that there are no bound states for a//3 < v — D /2 (for the special 
mode in the case a//3 = u — D/ 2 see below). For a//3 > v — D/2 a single bound state A = irj appears. 
The corresponding root r/ increases with increasing a//3 and for some critical value (3 = /3r = /3r(/c|^ ) 

one gets r] = k^j. Here, /3 r is the value of the Robin coefficient for which the root of the equation 

is equal to k^j. The stability of the vacuum state requires the condition 1//3 < 1//3r. The critical 
value for the Robin coefficient depends on the lengths of the compact dimensions, on the phases in 
periodicity conditions and on the mass of the field through the parameter u. Note that for a brane 
on AdS bulk with all dimensions being non-compact one has = 0 and all the modes with A = ir /, 
> 0, lead to the instability. Hence, in models with compact dimensions the stability condition, in 
general, is less restrictive. Assuming that a/ (3 > v — D/2, let us denote by u = u/ {a/ (3) the root of 
the equation () 2.291) . This root increases with increasing a//3 and does not depend on the location of 
the brane. The stability condition for the vacuum state is written as Uv\a/(3) < k^zo. From here it 
follows that, for fixed values of the other parameters, when the brane approaches the AdS boundary 
(zo decreases), started from the critical value zq = ul R \a/(3)/ k^j, the vacuum in the R-region becomes 
unstable. 

The coefficient in (12.261) is found from the normalization condition (12.201) making the replace¬ 
ment d(A — A') —» 5 rm i. By using the result for the integral involving the square of the Macdonald 
function [29] , one gets 


| C (b ),2 = (2tt) p a l D ril u (r)z 0 ) 

V q uj(rj)z 0 Kl(riz 0 ) 


(2.30) 


In deriving this expression we have used the relations 


K„(x) 


Bp 

Iu{x) ’ 


Bp(x 


2 + v 2 ) 



B 0 x 


K(x) 
K v {x) ’ 


(2.31) 


valid for x = rjzp, with r]zp being the solution of (12.281) . These relations are obtained by making use 
of (12.281) and the Wronskian relation for the modified Bessel functions. 

Now, in the presence of the bound state, the mode sum for the Hadamard function has two 
contributions. The first one comes from the modes with real A and is given by the expression (12.221) . 
The second contribution comes from the bound state. For the latter, by using the corresponding mode 
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functions (|2.26l) and the normalization coefficient from (12.301) . one finds 




2 a 1 d (zz') I 2 
(2ir) p V qZ Q 


E 


dk p e 


ik r Ax r 


X 


Iujyzp) cos(uj(y)At) 
K'uivz o) w(?7) 


K u (yz)K v (yz'). 


(2.32) 


The evaluation of the part coming from the modes with real A is similar to that we have described 
above. The difference arises in the step when we rotate the integration contour in (12.241) . Now, 
the integrand in this expression has poles A = ±ir/ on the imaginary axis, where y is the root of 
(|2.28p . After the rotation, the integration contour has to pass round these poles on the right by small 
semicircles. The integrals over the semicircles around A = iy and A = —iy give the residue at A = iy 
multiplied by 2iri. It can be seen that this residue term exactly cancels the contribution of the bound 
state in (12.321) . Hence, we conclude that the expression (12.251) for the Hadamard function is valid in 
the case of the presence of bound states as well. 

In addition to the modes discussed above, a mode may be present for which A = 0 and ui = k. 
For this mode the function Z v in (12.121) is a linear combination of z u and z ~ u . The part with z u is 
excluded by the normalizability condition and the mode functions have the form 



/2—v ^ik r x r ^ikt 


These modes are normalizable under the condition v > 1 and for the coefficient one finds 



/ i\ 2d—2 

[V - 1 ) Zp 

('.2ir) p V q a D ~ 1 k 


(2.33) 


(2.34) 


From the boundary condition on the brane it follows that the mode is allowed for the special value of 
the Robin coefficient determined from 


P/a = \/{y - D/2). 


(2.35) 


For v = D/2 this value corresponds to Neumann boundary condition and in this case the mode function 
does not depend on the coordinate z. An example of this special case is realized by a minimally coupled 
massless scalar field. Note that this special mode for a scalar field is the analog of the zero mode of 
the graviton in Randall-Sundrum 1-brane model [30] , 

For v > 1 and in the case of Robin boundary condition with (12.351) . the contribution of the special 
mode (12.331) should be separately added to the Hadamard function in formulas (12.22|) and (12.241) 
(but not to (12.251 ) . see below). Note that the mode function (|2.33l) can be written as A(pq CT ( x ') = 

where 

_ „D/2-v 

«(R)(*) = i)/ 2 , (2.36) 


and V§ ](T {x) are the mode functions for a massless scalar field in D-dimensional Minkowski spacetime 
with the spatial topology BP x T q . From here it follows that the contribution of the mode (12.331) to the 
Hadamard function is expressed as G(r\(x,x') = Q^(z)Q(^(z')G^ xTq (x, x'), where G^ xTq (x, x') 
is the corresponding function for a massless scalar held in D-dimensional Minkowski spacetime with 
the spatial topology R p x T q . By using the expression for the latter one gets 


^(R) {x, x') 


^(R)(z)Pl(R){z r ) ^ ik t Ax l K (p-i)/^(g) s p) 

(27t) P+1/2 Vq A ( Sp / fc(g) )(p-l)/2 ■ 


(2.37) 
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where the summation in the exponent goes over l 
defined s p = \/|Ax p | 2 — (At) 2 and 


p + 1,..., D — 1. In this expression we have 


k fq) = 2 k i = ( 2ym * + & d 2 / L i> ( 2 - 38 ) 

i=p +l l=p +1 

is the squared momentum in the compact subspace. 

It can be seen that the expression (12.251) for the Hadamard function is not changed by the pres¬ 
ence of the mode (12.331) . Indeed, under the condition (I2.35p . the contribution (12.371) of this mode 
is separately added to the right-hand sides of (12.221) and (12.241) . Now we should take into account 
that the integrand in (I2.24P with separate j has a simple pole at A = 0 and in the rotation of the 
integration contour this pole should be avoided by arcs of a circle of small radius. The contributions 
of the integrals over these arcs exactly cancel the contribution (12.371) of the special mode and, as a 
consequence of this, the representation (12.251) is not changed. 


2.3 Hadamard function in the L-region 

In the region between the brane and AdS boundary, y < yo, the mode functions still have the form 
(12.121) . From the normalizability condition it follows that for r^lwe should take 

Z v (Xz) = C a J u (Xz). (2.39) 


For the part of the solution with the Neumann function, the normalization integral diverges on the 
AdS boundary z = 0. In AdS/CFT correspondence normalizable and non-normalizable modes are 
dual to states and sources, respectively. The key feature of AdS spacetime is the presence of a timelike 
boundary at infinity where appropriate boundary conditions should be imposed in order to have well- 
defined dynamics 01- In the problem at hand, for the unique quantization procedure, in the range 
0 ^ v < 1, we need to specify the boundary condition on the AdS boundary. For a scalar field, the 
Dirichlet and Neumann boundary conditions are the most frequently used ones. The general class of 
allowed boundary conditions of the Robin type, has been discussed in |9j. Here, we choose (12.391) for 
all values of v ^ 0 that corresponds to Dirichlet condition on the AdS boundary in the case 0 ^ v < 1 
(note that the analytic continuation to the Euclidean section automatically selects this boundary 
condition EH). 

From the boundary condition on the brane at y = yo it follows that the eigenvalues of A are roots 
of the equation 

J»(Xz Q ) = 0, (2.40) 


where the notation with overbar is defined in accordance with (|2.18l) where now in (I2.19P 5 y = — 1. 
Hence, in the L-region the spectrum for A is discrete. First we consider the case when all the roots 
of (12.401) are real. Let us denote by x = 7 „, n = 1, 2,..., the positive zeros of the functions J u {x). 
Then, for the eigenvalues of A one has A = 'jn/zo (for a mode with purely imaginary A see below). The 
normalization coefficient is determined from (|2.20f) . with the replacement <5(A — A ; ) —>• 8 nn > and with 
the y-integration over (—oo,yo]. By using the standard result for the integral involving the square of 
the Bessel function [29j, we find 



\z 0 T u (\z 0 ) 
(2tt) P a D ~ 1 V q ZQL 0 ’ 


(2.41) 


with the notation 

T u (x ) = x[x 2 J' 2 {x) + {x 2 - u 2 )J 2 (x)}~ 1 . (2.42) 


Note that in the latter expression we could substitute xJ' u (x) = —AqJ u (x)/Bq. 
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Plugging the mode functions into the mode-sum (I2.11|l . for the Hadamard function we get the 
expression 


G(x, x') = 


2 a 


1 -D 


W) 


>\D /2 


(2vr yv q zl 


Y j dk p e iKAxr Y 


In 


cos (A Vyg/sff + k 2 ) 

V'yl/z o + A;2 


Tvip/n)JvisinZ/ " 0 ) Jv(SfnZ / 4 ))- 


(2.43) 


This expression involves the roots 7 n which are given implicitly and is not convenient for the further 
evaluation of the current density. In order to obtain more workable expression we apply to the series 
over n the summation formula [32] 


n= 1 


1 r°° 

2 Jo du ^ U ' ) 


27T Jo I v (u) 


x[e-»*/W + e m /N]. 


(2.44) 


The part in the Hadamard function coming from the first term in the right-hand side of (12.4411 gives 
the corresponding function in the geometry without the brane. As a result, we get the following 
decomposed representation 


G(x, x') 


r , ,, 4 0 ‘- D (zz'f 1 

Go{x, x ) - 


xA 


(2vrf +1 V q 
cosh(A t\Ju 2 — k 2 ) K v {uzo) 

Vu 2 - k 2 I v {uzo) 


Yj dk p e ikrAxr jT 
I u (uz)I v (uz'). 


du 


(2.45) 


Comparing with (12.25[1 . we see that the expressions for the brane-induced parts in the Hadamard 
function in the R- and L-regions are obtained from each other by the replacements I v (x) zz K u (x) 
and with the replacement /3 —» —/3 in the notations with overbars. Comparing with (12.4311 . we see 
the important advantages of the representation (12.4511 : (i) the contribution of the brane is manifestly 
separated, (ii) the explicit knowledge of the zeros 7 n is not required, (iii) the integrand in (12.451) is 
monotonic instead of the oscillatory behavior in (12.4311 and (iv) the representation (12.4511 holds in the 
presence of the bound state and of the special mode (see below) as well. In particular, the second and 
third points are important in the numerical evaluation of the vacuum currents. 

Depending on the value of the Robin coefficient /?, the equation (12.4011 can have purely imaginary 
roots A = ir /, > 0. In this case, for the stability of the vacuum state we should assume that rj < k^°j. 

From the boundary condition on the brane it follows that the allowed values for r/ are roots of the 
equation I u (r]Zo) = 0 which is written in the explicit form as 


(D/2 + u — a//3) I u (u ) + ul u+1 (u) = 0, 


(2.46) 


where u = rjzo- From here we conclude that the modes under consideration are absent in the case 
a/13 < D/2 + v (the special mode for the case a/j3 = D/2 + v is discussed below). For a//3 > D/2 + v, 
the equation (12.4611 has a single positive solution, u = iti L ^(a//3), which increases with increasing 
a/13. Started from the critical value of a/j3, denoted here by o//3l and determined from the condition 
ul L \a//3 l) = zq , the vacuum becomes unstable. Note that the critical values /3 r and /3 l of the 
coefficient in Robin boundary condition are different for the R- and L-region. As a result of this, 
there are values of /3 for which the vacuum is stable in the one region and unstable in the other. We 
see that, under the condition a//3 > D/2 + 17 when the location of the brane approaches the AdS 
boundary, started from the critical value zq = (a//3)/, the vacuum state becomes unstable. 
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The summation formula (|2.44l) is valid also in the presence of purely imaginary roots if we add to the 
left-hand side the contribution from the corresponding modes. This contribution has the form (12.431) 
with the replacement —» ir)Zo and omitting the summation over n. As a result, the representation 
(|2.45l) is valid in the presence of purely imaginary roots as well. 

Similar to the case of the R-region, under the condition 

P/a = 1/ (D/2 + v) , (2.47) 


there is a A = 0 mode with the mode function 


= C {L )Z 


D/2+v ik r x r ^ikt 


where 


Cfr , = 


(y + l)a 


l—D 


(L) ( 27 t) p V qZ Q U+2 k' 


(2.48) 


(2.49) 


The contribution of this mode to the Hadamard function is obtained by taking into account that 

= %)(^(m) ct ( x )’ with 


V/2+v 

n (L) {z) = y/2(v + l)a (1 ~ £>)/2 u+1 , (2.50) 

z o 

and is obtained from (12.371) by the replacement J4( R )(z) —» f2( L )(z). This contribution should be added 
to the right-hand side of (12.43[) . The representation ()2.45l) is not changed. 


3 Vacuum currents 


By using the expressions for the Hadamard function, from formula (12.101) we can see that the VEVs 
of the charge density and of the components of the current density along non-compact dimensions 
vanish: 

(j l ) =0, 1 = 0,1,..., p. (3.1) 

Of course, the latter property for the spatial components is a direct consequence of the problem 
symmetry under the reflections x l —>• — x l . For the component along the Zth compact dimension one 
finds the decomposition 

(j l ) = (j l )o + (j l ) b , l=p+l,...,D-l, (3.2) 

where (j l )o is the corresponding VEV in the absence of the brane and the part (j l )b is induced by the 
brane. 

The contribution (j l )o is investigated in [19] and is given by the expression 


(j l 


'o = 


lZ)(D+m ^ n i sh <ami) Y, cos(^cw) 


+XX l * 2 /( 22;2 ))’ 

i 

where n ? _i = (n p+u ... ,n;_i,rq + i,.. .,n D -i), and 


_ e-WQ^x) 
(x 2 — l)^/ 2 ’ 


(3.3) 


(3.4) 
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with Q&(x) being the the associated Legendre function of the second kind. Near the AdS boundary, 
z —>• 0, the current density (13.31) behaves as z D+2u+2 and near the horizon the leading term in the 
asymptotic expansion is given by (j l )o ~ ( z / a ) D+1 {j l )^>+ i xTg - Here, 


(M M ) 

/ RP+l XT’ 


2eL, 


r((Z3 + l)/ 2 ) 

7 t (- D + 1 )/ 2 


oo 

m sm(aini) ^ 

n ; =l n 5 _i 


COs(£ OCiUi) 


(3.5) 


is the VEV of the current density for a massless scalar field in Minkowski spacetime with spatial 
topology _R P+1 x T q . 

Similar to the case of the Hadamard function, we shall consider the brane-induced contribution in 
the YE Vs of the current density for the R- and L-regions separately. 


3.1 R-region 

In the R-region the brane-induced contribution in (13.21) is obtained from the corresponding part of the 
Hadamard function in (12.251) . By using the relation 



uf{u)du 

Vu 2 — k 2 


( P / 2 ) 

2 r((p + 1 )/ 2 ) 


duu p f(^Ju 2 + kf q) ), 


we find the expression 


(A = 


eC p z D+2 
2 P~ l a D + l V q 




Iy(z 0 x) 

K v {z 0 x) 


Kl(zx), 


(3.6) 


(3.7) 


with the notation 

Tr-(P+ 1)/2 

C p = —-——. (3.8) 

p r((p + 1 )/ 2 ) 1 ; 

From (13.71) we see that the brane-induced contribution to the current density along the Zth compact 
dimension is an odd periodic function of the phase oy with the period 2ir and an even periodic function 
of the remaining phases dj, i ^ l, with the same period. By taking into account the relation (12.71) . 
we conclude that the VEV of the current density is a periodic function of the magnetic flux with the 
period equal to the flux quantum. 

The charge flux through the (D— l)-dimensional spatial hypersurface x l = const, having the normal 
ni = a/z, is given by the quantity ni(j l ). From ()3.7I) it follows that the corresponding brane-induced 
contribution, ni(j l )b , depends on the lengths of compact dimensions and on the coordinate 2 in the 
form of the ratios Li /zq and z/zq. The latter is expressed in terms of the proper distance from the 
brane, y — yo, as z/z$ = e^ y ~ yo ^ a . 

Let us first consider the flat spacetime limit of (|3.7I) . corresponding to the limiting transition 
a —>• oo for fixed values of y and yo- I n this limit, the order of the modified Bessel functions in (13.71) 
is large. Changing the integration variable to i -> vx, we use the corresponding uniform asymptotic 
expansions (see, for instance, m)- By taking into account that in the limit under consideration 
z ~ a + y and zq « a + yo, to the leading order we get 


(A 



*&) - m3) V e -* ( v-»>!g±i 


(3.9) 


The expression in the right-hand side coincides with the boundary-induced part of the current density 
for the geometry of a single Robin plate in (D + l)-dimensional Minkowski spacetime with spatial 
topology R p+1 x T q (see [23]). 
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For a conformally coupled massless field one has v = 1/2 and the modified Bessel functions in (13.71) 
are expressed in terms of the elementary functions. In this case the expression for the total current 
density takes the form 


(j l ) = ( z/a) D+1 


(ji\( M ) 

\J /rp+IxTI 


+ y kl r dx (* 2 - m£±I 

4.) <a) -si*- 1 


(3.10) 


where (the notation with the lower sign is employed below) 


Pm = 


Pz 0 /a 


l±(D-l)/3/(2ay 


(3.11) 


The right-hand side of (I3.10j) . divided by the conformal factor ( z/a) D+1 , coincides with the current 
density in the corresponding problem on Minkowski bulk with the plate at z = zo (compare with 
(13.91) 1 on which the field obeys the Robin boundary condition (12.81) with the replacement (3 /3^. 

The difference of the Robin coefficients in the two conformally coupled problems is related to the fact 
that this coefficient is not conformally invariant. 

At large distances from the brane compared with the AdS curvature radius, y — yo ^$> a, one has 
z 2> zq. In addition, assuming that z^> Li, we can see that the dominant contribution to the integral 
in (13.71) comes from the region near the lower limit and the contribution of the mode with a given n 5 
is suppressed by the factor e~ 2zk ^. Under the condition \oti\ < 7r, assuming that all the lengths L,- 
are of the same order, the main contribution comes from the term with m = 0, i = p + 1,... ,D — 1, 
and to the leading order we find 


ez D (p 1)/2 d^ ( ( ° )(p 1)/2 Uzok^) _ 2zk m 
2P+^(p-m a D + i VqLl RAzok [°J) 6 


(3.12) 


This asymptotic corresponds to points near the AdS horizon. As we have already mentioned, in this 

limit, for the boundary-free part one has (j l )o ~ ( z /a) D+1 (j l )^\ lxTq , where {j l )^+i xTq is the current 

density for a massless scalar field in (D + l)-dimensional Minkowski spacetime with spatial topology 
RP +i x T q ( see 

(13.51) 1 and with the lengths of the compact dimensions Li, i = p + 1,..., D — 1. From 
here we conclude that near the horizon the boundary-free part dominates in the total VEV. 

For fixed values of z and Lj, when the location of the brane tends to the AdS boundary, zq —>• 0, 
to the leading order, from (13.71) one finds 


/ -i\ ^ 4eC p z D+ ~Z q v Aq + vBq 2 v+p+i 

U )b - 2 2u+ PvT 2 (v)a D+1 V q A 0 - uB 0 1 ^ 

/ OO 

dxx 2u+1 (x 2 - l )(P~ l)/2 K 2 (zk {q) x), (3.13) 

and the VEV vanishes as z$ v . 

Now, let us consider the limit when the length of the Zth dimension is much smaller than the 
lengths of the other compact dimensions, Li <C Li. In this case, in (13.7p the dominant contribution to 
the sum over m, i = p + 1,... ,D — 1, i ^ l, comes from large values of \rii\ and we can replace the 
summation by the integration in accordance with 


J2 /( fc u-i)) 

n q -i 


2 2- % -( ? -1)/2 V t 

L,r(( g -i)/2) 



duu q 2 f{u), 


(3.14) 


where k 2 q _^ = k 2 q ^ — k 2 . By making this replacement 
integration variable w according to x = Jw 2 + u 2 + k 2 . 


in (13.71) . instead of x we introduce a new 
Then, introducing polar coordinates in the 
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plane ( u,w ), the integral over the angular variable is expressed in terms of the gamma function. As 
a result, to the leading order we get 


(j l 


eC D . 2 z D + 2 
2 D - 3 a D+ 1 L, 


+oo 


^ ki / dxx(x 2 — kf) 


2 , 0=3 I u (z 0 x) 


ni =—oo 


%i 


K u (z 0 x) 


I<l(zx). 


(3.15) 


The expression in the right-hand side coincides with the brane-induced contribution in the model with 
a single compact dimension of the length Li (q = 1, p = D — 2). 

If in addition to L; <C Li one has L/ zo, the arguments zqx of the modified Bessel functions in 
(13.151) are large. By using the corresponding asymptotic expressions [33], after the integration over x 
we find 


tA 


(1 - 25 0Bo )e(z/a) D+1 
2 D ~ 2 7T D / 2 Li ( Z — Zq ) D / 2 1 


+00 

Y ki\ki\ D/ 2 ~ 1 K D/2 _ l {2 (.z - zo) \k t \). 

m=—oo 


(3.16) 


Here, for non-Dirichlet boundary conditions we have assumed that |/3|/a Li/zq. From (|3.16|) it 
follows that the brane-induced contribution is located near the brane within the region z — zq < Li 
and has opposite signs for Dirichlet and non-Dirichlet boundary conditions. At distances z — zq Li 
it is suppressed by the factor e~ 2 ( z - z o) a i/ L i . Note that, in the limit Li <C Li,z, for the boundary-free 
part one has the asymptotic 


(j‘)o 


2er((g+ l)/2) 

n ( D + 1 )/ 2 {a/z) D+l L^ 


E 

n= 1 


sin {pun) 


n 


D 


(3.17) 


The expression (13.71) is not convenient for the investigation of the current density behavior near 
the brane. To this aim, a more convenient expression for the VEV of the current density is obtained 
by using the representation (IA.3D for the Hadamard function. After the integrations over w and k p , 
we get the following result 


( 3 l ) 


4 ea D 1 z D+2 ^ sin {non) f°° ,, . 


g 2 (\z 0 :Xz) / 2 - 

J^(A= 0 ) + F„=(A^ 0 ) a ' ,/2+l( 'V 


(3.18) 


where we have defined the function 

g v {x) = x u K v {x). (3.19) 

In the presence of a bound state, its contribution should be added to (|3.18D separately. In the model 
with a single compact dimension x l of the length L and with the phase a, from ([3.181) one finds 


4ea D l z D+2 ^ sin (na) f°° 

(2n)D/2 ld-i A, J o 

gl(Xzo 1 Xz)g D / 2 {nLX) 
J 2 (Xz 0 )+ Y 2 (Xz 0 ) 


(3.20) 


where we have substituted p = D — 2. 

An important result which follows from (13.181) is that the VEV of the current density is finite on 
the brane. The corresponding value is directly obtained from (13.181) putting z = zq and by taking 
into account that g u (u,u) = 2 Bq/it. For Dirichlet boundary condition both the current density and 
its normal derivative vanish on the brane. The finiteness of the current density is in clear contrast to 
the behavior of the VEVs for the field squared and the energy-momentum tensor which suffer surface 
divergences. For example, the VEV of the field squared diverges as 1 /{z — zq) D 1 • The feature that 
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the VEV of the current density is finite on the brane could be argued on the base of general arguments. 
In quantum field theory the ultraviolet divergences in the VEVs of physical observables bilinear in 
the field are determined by the local geometrical characteristics of the bulk and boundary. On the 
background of standard AdS geometry with non-compact dimensions the VEV of the current density 
in the problem under consideration vanishes by the symmetry. The compactification of the part of 
spatial dimensions to g-dimensional torus does not change the local bulk and boundary geometries 
and, consequently, does not add new divergences to the expectation values compared with the case of 
trivial topology. 

Under the condition (12.351) . the contribution of the special mode (12..' >3 1) has to be added to the 
right-hand side of (13.181) (representation (13.71) is not changed). This contribution is obtained from 
the corresponding part (|2.37l) in the Hadamard function and is related to the current density for a 
massless scalar field in D-dimensional Minkowski spacetime with the spatial topology RP x T q by the 
formula (j l )(R) = ^(r )( z )(j l )ftP xT < n where the factor is given by the expression (12.361) . By 

using the expression for OOrpxT? from [20], one gets 


(j l ) (R) 


8e (v — 1) Zq U 2 z d 2v ^ sin ( noq ) 
(27r)(P+ 3 )/ 2 a D - 1 U Q Lf ^ n?+ l 


E flp/2+l(™fr^(g-l))- 

n q ~ i 


(3.21) 


Recall that, as the necessary condition for the presence of this contribution we have v > 1. In the 
model with a single compact dimension this gives 


(j l ){ R) 


4e {u - 1) T(D/2) 2 „_ 2 D _ 2v ^ sin (non) 
TT D / 2 a D ~ 1 L D ~ 1 Z ° Z ^ n D ~ l 

n —1 


(3.22) 


Note that in the case v = D/2, the current density from the special mode does not depend on z. 

Let us consider the behavior of the current density in the limit when both the location of the brane 
and the point of observation are close to the AdS boundary, zo,z <C L*. Under this condition, the 
arguments of the Bessel functions J v (x) and Y v (x ) in ()3.181) are small and we use the corresponding 
asymptotic expressions. To the leading order, the integral over A is expressed in terms of the Macdonald 
function and we get the expression 

„ ea~ D ~ l z D ~‘ 2u+2 L ~^ p ~ 2v ~ 2 ( 2v Aq + vB 0 2 A 2 

> ~ 2P/2+^-l 7r p/2+i r ( i/ + 1 )V q V Ao - ) 

x ES^SrEw2«+2("i.Vi))' (3-23) 

ti=l n g _i 

The part with the term z 2v in the brackets gives the asymptotic for the VEV in the absence of the 
brane. 

The representation (13.181) is also well-suited for the investigation of the asymptotic behavior in the 
limit of large L; compared with the other length scales of the model. In this limit the argument of the 
function g tl (x) in (13.181) is large and we can use the asymptotic expression <y /t (x) ~ /2x fJ ’~ 1 / 2 e~ x , 

for i> 1. The dominant contribution to the integral comes from the region near the lower limit. Let 
us denote by the lowest value of k( q _ i), = min(fc( g _ 1 )). For |dj| < tt this value is realized 

by the mode with n* = 0, i ^ l, and we have 

4°-i)= E ^/ L i- ( 3 - 24 ) 

i=p+l,jt=l 
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Two cases should be considered separately. For ^ 0, in the series over n the leading contribution 

comes from the n = 1 term and we get 


(?) 


2 ea~ D ~ 1 z D ~ 2l,+2 k^^ 3 ‘^ 2+u 
7i-(P+i)/ 2 r(i/ + l)V g (2L i )(P+ 1 )/2+^ 
( 2v A 0 + vB 0 ~ 2v \ 2 sind; 

V Ao-„Bo~°) 


(3.25) 


In this case the current density is exponentially small. Note that (13.251) is also obtained from (13.231) in 
the limit of large L;. For = 0, by using the standard integral involving the Macdonald function 

[29] , for the leading term in the asymptotic expansion one finds the expression 


4eT (p/2 + v + 2) a D l z D 2v+2 

nP/ 2 +iT(i' + 1 )V q L/ +2u+2 

( 2v A 0 + vBo 2u\ 2 V'' sin 

X V~ Aq- vB 0 Z ° J ^ n p+2u+3 ' 


(3.26) 


Now the decay of the current density with increasing Li is power law for both massless and massive 
fields. This result for massive fields is in contrast to the corresponding behavior of the current density 
in Minkowski bulk. In the latter geometry the decay of the current density is exponential with the 
factor e~ mLl . 

For points outside the brane, z > zo, another expression for the current density is obtained from 
the representation (IA.5I) . After evaluating the integrals, one finds the following result 


O') 


, u 4ea 1 D L l p z d+2 Z4 sin (ncq) ^ 
0 >0-TFZw^ZTT?- 2^ „P +1 


( 2ir)P/ 2 + 1 V i 


nP + 1 

9 71=1 n,j_l ^ K (q-1) 


dx 


XX r!(xZ 0 \ K v( XZ ) W p/ 2 +l( nL l\J x2 - k fq- !))• 


(3.27) 


with the function 

w v (x ) = x v J v {x). (3.28) 

In the absence of the bound states, the equivalence of the representations (|3.7I) and (13.271) for the 
brane-induced contribution is directly seen by using the formula 


+°° 2X/ °° /*oo 

MN) = — ^ sin(nd() / dx xsm(nLix)g(x). (3.29) 

n;=—oo n =1 ' 

The latter relation follows from the Poisson’s resummation formula (see also {23]). 

In what follows, all the graphs are plotted in the D = 4 model with a single compact dimension of 
the coordinate length L and with the phase d, for a minimally coupled (£ = 0) massless scalar field. 
For the corresponding value of the parameter v one has v = D/ 2 = 2. 

On the left panel of figure [2] we have depicted the current density as a function of the phase a 
for fixed values of the parameters zq/L = 1, z/zq = 1.2. The graphs are plotted for Dirichlet (D), 
Neumann (N) and for Robin (with /3/a = —1, the number near the curve) boundary conditions. The 
dashed curve presents the current density in the same model when the brane is absent. As is seen, 
depending on the boundary condition, the presence of the brane leads to the increase or decrease of 
the current density. Note that for the example considered one has v — D/ 2 = 0 and, in accordance 
with ()2.35l) . for Neumann boundary condition there is a special mode (12.331) the contribution of which 
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Figure 2: The VEV of the current density as a function of the phase in the periodicity condition (left 
panel) for D = 4 AdS space with a single compact dimension and for Dirichlet, Neumann and Robin 
(/3/a = —1) boundary conditions. The graphs are plotted for zq/L = 1, z/zq = 1.2. The right panel 
displays the ratio ni(j l )/(j l ) m as a function of the coefficient in Robin boundary condition for fixed 
values a = 7t/2, z/zq = 1.2, zq/L = 1. 


is given by (13.211) . This contribution should be added to (13.201) . Alternatively, in the numerical 
evaluation we may use the representation (|3.7I) which holds in the presence of the special mode as 
well. We have numerically checked that both these ways of the evaluation give the same result. 

The right panel of figure [2] presents the ratio ni(j l )/(j l ) m as a function of the coefficient in Robin 
boundary condition, measured in units of AdS curvature radius. Here, 


(j 1 )m 


2eT((D + l)/2) 
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sin (cm ) 


n 


D 


(3.30) 


is the current density for a massless scalar field in (D + l)-dimensional Minkowski spacetime with 
topology it! 0 ” 1 x S' 1 and with the length of the compact dimension aL/z. Note that the latter is 
the proper length of the compact dimension in AdS spacetime measured by an observer with a given 
z. The graph is plotted for fixed values a = n/2, z/zq = 1.2, zq/L = 1. The vertical dotted curve 
corresponds to the critical value /3r/o ~ 1.31. In the region 0 < /3 < /3r the vacuum is unstable. The 
horizontal dashed curves correspond to Dirichlet and Neumann boundary conditions. As expected, 
in the limits j3 — >• 0 and (3 —>• oo the results for Dirichlet and Neumann conditions are obtained. For 
Dirichlet boundary condition the current density takes its minimal value (the minimal absolute value 
for negative a). 

In figure [3] we have plotted the ratio ni(j l )/{j l )m as a function of z/zq for Dirichlet (left panel) 
and Neumann (right panel) boundary conditions on the brane. The numbers near the curves are the 
values of zq/L. For the phase we have taken the value a = ir/2. As it has been already mentioned, for 
Dirichlet boundary condition the current density vanishes on the brane. From the asymptotic analysis 
given above it follows that on the horizon, z/zq —» oo, one has ni(j l )/(j l )yi —» 1. This behavior is seen 
from the graphs. 

Figure |4] displays the dependence of ni(j l )/(j l ) m on z/zq in the case of Robin boundary condition 
for several values of (3/a (numbers near the curves). For (3/a = —2/(D — 1) the normal derivative 
d z (ni(j l )/(j l ) m) vanishes on the brane. This can also be directly seen from the analytic expression 
by taking into account that vd v g v (u,v)\ v = u = —2Ao/ir. 

From the results derived in this section we can obtain the current density in /^-symmetric 
braneworld models of the Randall-Sundrum type with a single brane. In the original Randall-Sundrum 
1-brane model [30] the universe is realized as a ^-symmetric positive tension brane in 5-dinrensional 
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Figure 3: The dependence of the quantity ni(j l ) / (J 1 )m on z/zq for Dirichlet (left panel) and Neumann 
(right panel) boundary conditions. The graphs are plotted for fixed values of zq/L (figures near the 
curves) and for a = 7t/2. 



Figure 4: The same as in figure [3] for Robin boundary condition. The graphs are plotted for Zq/L = 1 
and the numbers near the curves correspond to the values of (3/a. 
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AdS spacetime. In this simplest variant the only contribution to the curvature comes from the negative 
cosmological constant in the bulk. However, most scenarios motivated from string theories predict 
the presence of other bulk fields, such as scalar fields. In addition, they predict also small compact 
dimensions originating from 10D string backgrounds. In a generalized (D + l)-dimensional version of 
the Randall Sundrum 1-brane model the line element is given by (|2.2I) with the warp factor e~^ y ~ y °^ a 
where yo is the location of the brane. The background geometry contains two patches y > yo of 
the AdS glued by the brane and related by the ^-symmetry identification y — yo <—> Vo ~ U- The 
corresponding spatial geometry in the case D = 2, embedded into a 3D Euclidean space is depicted 
in figure [5j 



Figure 5: D = 2 spatial geometry corresponding to the Randall-Sundrum 1-brane model with a 
compactified dimension. 

Because of the absolute value sign in the exponent of the warp factor, the Ricci scalar contains a 
contribution located on the brane, 

R = 4D5{y-y 0 )/a- D{D + l)/a 2 . (3.31) 

For non-minimally coupled scalar fields, this leads to delta-type terms in the field equation. An 
additional delta-type term may come from the boundary action of a scalar field of the form Sb = 
c f d D xdy \/\g\(p 2 5(y — yo), where c is the so-called brane mass term. The boundary condition for the 
mode functions is obtained by integrating the field equation near the brane. In a way similar to that 
used in puna mg, it can be seen that for fields even under the reflection with respect to the brane 
(untwisted scalar field) the boundary condition is of the Robin type with 

P = -l/(c + 2D£/a). (3.32) 

In particular, for a minimally coupled field with the zero brane mass term the boundary condition 
is the Neumann one. For fields odd with respect to the reflection (twisted fields) the boundary 
condition is reduced to the Dirichlet one. Now, in the ^-symmetric model the integration over y 
in the normalization integral (12.201) goes over the interval (—oo,+oo). As a result the square of the 
normalization coefficient contains an additional factor 1/2 compared to the one we have obtained for 
the R-region. Hence, the expressions for the VEV of the current density in the generalized Randall- 
Sundrum 1-brane model with compact dimensions are obtained from those given above in this section 
with an additional factor 1/2 and with the Robin coefficient (13.321) for untwisted fields and with 8 = 0 
for twisted fields. 


3.2 L-region 


Now we turn to the current density in the L-region. By using the expression (12.451) for the Hadamard 
function, the current density in this region is decomposed as (13.21) with the brane-induced part 
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(3.33) 


20 













Under the condition r] < k !^\, this representation is valid in the presence of the mode with A = ir /, 
where rj is the zero of the function I u (zor)), and also in the presence of the special mode (12.48[) . For 
large values of AdS radius a, in a way similar to that for the R-region, we can see the limiting transition 
of the expression (13.331) to the corresponding formula for a plate in Minkowski bulk. 

For a conformally coupled massless held, the expression of the total current density takes the form 
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(3.34) 


with dehned by (13.111) . Here, the first term in the figure braces and the part with the first term in 
the square brackets come from (j l )o- The expression on the right of (13.341) . divided by the conformal 
factor ( z/a ) D+1 , coincides with the current density in the region between two plates on Minkowski 
bulk with Dirichlet boundary condition on the left plate and Robin condition (12.81) . with (5 —>• /3A, 
on the right one (see [24] for the problem with Robin boundary conditions on both plates). The fact 
that the problem with a single brane in AdS bulk in the L-region is conformaly related to the problem 
with two plates in Minkowski bulk is a consequence of the boundary condition we have imposed on 
the AdS boundary. 

The asymptotic behavior of the VEV near the AdS boundary, z —>• 0, is directly obtained from 
(13.331) . by using the expression of the modified Bessel function for small arguments. To the leading 
order we get 




r)l-2i/—p /~i „D+2u+2 _ 

2 y eO p Z _, r2u+p+l 

a D+1 V q T 2 (v + 1) jp w) 


/ OO 

dx x 2v+1 {x 2 


1 y P -l)/2 Ku ( Z °k(q) X ) 
I u (z 0 k [q) x) ’ 


(3.35) 


and the brane-induced contribution vanishes as z D+2u+2 . Recall that near the AdS boundary the part 
(j l )o in the VEV of the current density behaves in a similar way and, hence, on the AdS boundary 
the ratio of the brane-induced and boundary-free contributions tend to a finite limiting value. 

In the limit when the brane tends to the AdS horizon, zq —» oo, the argument zqx of the modified 
Bessel functions is large. By using the corresponding asymptotics and by taking into account that the 
dominant contribution to the integral in (13.331) comes from the region near the lower limit, we can see 
that the leading order term is given by 
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(3.36) 


Hence, for a fixed value of z, when the brane location tends to the AdS horizon, the brane-induced 
contribution is exponentially suppressed. 

If the length of the It h dimension is much smaller than the lengths of the remaining compact 
dimensions, Li <C Li, in a way similar to that for the R-region, we can see that, to the leading order, 
the brane-induced contribution coincides with the corresponding quantity in the model with a single 
compact dimension of the length L[. The expression for the latter is obtained from the right-hand 
side of (|3.15|) by the replacements I v K v . If in addition L[ <C z, the corresponding asymptotic 
expression is given by the right-hand side of (13.161) with zq — z instead of z — zq, and the brane-induced 
contribution is concentrated near the brane in the region zq — z < L[. 
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The representation (13.331) is not well adapted for the investigation of the asymptotic near the brane. 
A more suitable representation is obtained by using the formula (1A.6I) for the Hadamard function: 
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(3.37) 


Unlike to the representation (13.331) . in the presence of a bound state, its contribution must be addi¬ 
tionally added to (13.371) . The latter is obtained from the right-hand side of ((3.37)1 by the replacement 
7 i —>• ir/ZQ and omitting the summation over i. The corresponding representation is valid under the 
condition g < From (13.371) we conclude that the VEV of the current density is finite on the 

brane. Similar to the VEV in the R-region, the current density and its normal derivative vanish on 
the brane for Dirichlet boundary condition. Another representation is obtained from (13.271) by the 
replacements I u ^ K u and with the replacement (5 —>• —/3 in the notations with overbars.. 

For the value of the Robin coefficient ()2.47l) there is a special mode (|2.48l) . As we have seen above, 
the corresponding contribution to the Hadamard function is expressed in terms of the Hadamard 
function for a massless scalar field in D-dimensional Minkowski spacetime with the spatial topology 
R p x T q . By using this relation, for the contribution of the special mode to the current density in the 
L-region we find 
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In the case of a single compact dimension this simplifies to 
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(3.39) 


Under the condition (12.471) . the contribution (13.381) should be added to the right-hand side of (13.371) . 

For the investigation of the asymptotic behavior for the contribution of the modes with A = 7 n /a 
in the limit of large values of L;, compared with the other length scales, it is convenient to use the 
representation (13.371) . By using the asymptotic expression of the function g^(x) for large values of the 
argument, we can see that the dominant contribution comes from the lowest mode with i = 1 and to 
the leading order one finds 
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with x = Li \J~ii /vj + ^’(V/- 2 i)' Hence, unlike to the R-region, the decay of the current density is 
exponential for both cases = 0 and 7 ^ 0. Under the condition (|2.47l) one has the additional 

contribution (13.381) from the mode with A = 0. For 7 ^ 0, in the limit of large Li, this contribution 

decays exponentially, as e Llk (i~ 1 ). In the case = 0 the decay is power law, like 1/L^ +1 . In both 

cases the contribution of the special mode dominates in the total VEV. 

Now let us consider the asymptotic of the current density (13.371) when the brane is close to AdS 
boundary, zq <C L t . The dominant contribution comes from large values of \rij\. j / l, and we can 
replace the summation over n„_i by the integration in accordance with (13.141) . After the integration 
over we get 
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The right-hand side presents the brane-induced contribution in the model with a single compact 
dimension of the length L/. In the limit under consideration the argument of the function g D / 2 (x) is 
large. By using the corresponding asymptotic expression, we see that the main contribution comes 
from the term with n = 1 , i = 1 with the result 
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(3.42) 


For the Robin boundary condition with (12.471) . the asymptotic for the contribution from the special 
mode is directly obtained from (13.381) . This contribution behaves as {zq/Li) d ~ 2 {z/ zq) d+2u and the 
corresponding decay, as a function of zq (for a fixed z/zq) is power law. 
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Figure 6: The current density in the L-region as a function of d for Dirichlet, Neumann and Robin 
(with /3/a = —1/2) boundary conditions (left panel) and as a function of the coefficient in Robin 
boundary condition (right panel). The graphs are plotted for zq/L = 1, z/zq = 0.8. On the left panel 
the dashed curve corresponds to the current density in the geometry without the brane. The vertical 
dotted line on the right panel corresponds to the critical value /3 c /a ~ 0.228. 

The left panel of figure[6]displays the vacuum current density in the L-region for Dirichlet, Neumann 
and Robin (with /3/a = —1/2) boundary conditions as a function of the phase in the quasiperiodicity 
condition along the compact dimensions. Recall that in the numerical evaluations we consider aD = 4 
minimally coupled massless field. In this case v = 2 and the modes with purely imaginary A and the 
special mode with A = 0 are absent for a/ f3 <4. The graphs are plotted for zq/L = 1, z/ zq = 0.8. The 
dashed curve corresponds to the current density in the absence of the brane. The right panel of figure 
[6] presents the dependence of the ratio ni(j l ) / (j l ) m on /3/a. The horizontal dashed lines correspond 
to Dirichlet and Neumann boundary conditions. The vertical dotted line corresponds to the critical 
value /3 c /a ~ 0.228. In the region 0 < /3 < (3 C the vacuum is unstable. 

In figured for a = 7r/2, we show the dependence of the ratio ni(j l )/(j l ) m on z/zq for Dirichlet (left 
panel) and Neumann (right panel) boundary conditions for separate values of zq/L (numbers near the 
curves). For Dirichlet condition the current density vanishes on the AdS boundary and on the brane. 
As it follows from the asymptotic (13.351) . near the AdS boundary the charge flux density ni(j l ) behaves 
as z D+2u+1 . For the Minkowskian VEV with the length of the compact dimension aL/z, equal to the 
proper length on the AdS bulk, one has {j l )m oc z d+1 . Hence, the ratio plotted in figure [7] vanishes 
on the AdS boundary as z 2v . Similar graphs in the case of Robin boundary condition are presented 
in figure [8] for several values of /3/a (numbers near the curves). Note that for all the examples in the 
L-region there are no modes with purely imaginary A. 
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Figure 7: The ratio ni(j l )/(j l ) m versus z/zq for Dirichlet (left panel) and Neumann (right panel) 
boundary conditions. The figures near the curves correspond to the values of zq/L and for the phase 
we have taken ot = 7t/2. 



Figure 8: The same as in figure[6]for the fixed value z/zq = 1 in the case of Robin boundary condition. 
The numbers near the curves correspond to the values of /3/a. 
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4 Conclusion 


We have studied the effects induced by a brane, parallel to the AdS boundary, on the VEV of the 
current density for a massive charged scalar field with an arbitrary curvature coupling parameter. 
The background geometry under consideration is obtained from the (D + l)-dimensional AdS one, 
described in Poincare coordinates, by a toroidal compactification of a part of spatial dimensions. 
Along compact dimensions the field operator obeys quasiperiodicity conditions with general constant 
phases. We have also assumed the presence of a constant gauge field. By the gauge transformation the 
problem is reduced to the one with a zero gauge field. In the new gauge the phases in the periodicity 
conditions are shifted by an amount determined by the ratio of the magnetic flux, enclosed by a 
compact dimension, to the flux quantum. On the brane and on the AdS boundary the field operator is 
constrained by Robin and Dirichlet boundary conditions, respectively. We consider a non-interacting 
quantum field and the Hadamard two-point function contains all the information about the properties 
of the vacuum state. Though the background geometry is homogeneous, the extrinsic curvature tensor 
of the brane is nonzero and, as a consequence of this, the regions on the right (R-region) and on the 
left (L-region) of the brane are not physically equivalent. 

In the R-region the spectrum of the quantum number A is continuous and the mode functions 
obeying the boundary condition on the brane have the form (12.1211 with the radial function (12.1611 . In 
addition to these modes, under the condition a//3 > v—D/ 2, one can have a bound state with the mode 
function (12.2611 and with the eigenvalue for r] being the root of the equation (|2.28p . The root should 
be constrained by r] < k^j, where k^j is defined as (12.2711 . For the modes with r/ > k^j the energy 
becomes imaginary and they lead to the instability of the Poincare vacuum state. In addition to the 
parameters v and the stability of the vacuum depends on the location of the brane. Depending 
on the value of the Robin coefficient, the approaching of the brane to AdS boundary can lead to the 
instability. For the value of the Robin coefficient given by (12.3511 there is a special mode (12.3311 for 
which A = 0. For a minimally coupled massless scalar field the special value of the Robin coefficient 
corresponds to Neumann boundary condition and the mode function does not depend on the radial 
coordinate z. This is an analog of the graviton zero mode in Randall-Sundrum 1-brane model. The 
Hadamard function in the R-region is given by the expression (12.251) with the second term in the 
right-hand side being the brane-induced contribution. We have shown that this expression holds in 
the presence of bound states as well. Alternative representations for the Hadamard function in the 
R-region, (lA.3j) and (1A.5I) . are derived in Appendix, by using the summation formula (1A.1I1 . 

In the L-region the spectrum for A is discrete and its eigenvalues are roots of the equation (12.4011 . 
The corresponding expression for the Hadamard function contains series over these roots. Another 
expression, in which the explicit knowledge of the eigenvalues for A is not required is obtained by 
making use of the generalized Abel-Plana formula (|2.44l) . This allows us to extract manifestly the 
brane-induced contribution, given by the second term in the right-hand side of (I2.45p . For a//3 > 
D /2 + v there is a mode with purely imaginary A = ir\ which is the root of the equation (|2.46l) with 
x = rjZQ. For the stability of the vacuum state one needs to have r] < The stability condition 
depends also on the position of the brane. The shift of the brane to the direction of the AdS boundary 
gives rise to the instability in the initially stable vacuum state. An alternative expression for the 
Hadamard function in the L-region is given by (j A.6I) . 

In both the R- and L-regions the VEV of the current density is decomposed into the boundary-free 
and brane-induced contributions. For both these contributions, the component of the current density 
along the Zth compact dimension is an odd periodic function of the phase cq and an even periodic 
function of the phases cij, i / l, with the period equal to 27 t. In the R-region the brane-induced 
contribution is given by (13.711 . We have checked that for large values of the AdS curvature radius the 
leading term in the corresponding asymptotic expansion coincides with the boundary-induced part 
of the current density in Minkowski spacetime with topology RP +l x T q in the presence of a single 
Robin plate. In the case of a conformally coupled massless field the current density coincides, up to 
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the conformal factor (z/a) D+1 , with the corresponding quantity in Minkowski bulk with the plate at 
z = zq on which the field obeys the Robin boundary condition with the coefficient given by (13.111) . 
For points near the AdS horizon (z is large compared with the other length scales), the brane-induced 

—2 zk^ 

contribution is suppressed by the factor e (<?). In the same limit the boundary-free part behaves 
as ( z/a) D+l and it dominates in the total VEV. In the limit when the brane approaches to the AdS 
boundary, for fixed values of z and L*, the brane-induced contribution tends to zero like Zq". For 
the investigation of the near brane asymptotic of the vacuum current it is more convenient to use the 
representation (13.181) . In the presence of a bound state its contribution should be added separately 
to the right-hand side. For Robin boundary condition with (12.351) an additional contribution, given 
by ()3.21|) . comes from the special mode (12.331) . An important conclusion which follows from the 
representation (13.181) is that the current density is finite on the brane. This behavior is in clear 
contrast with that for the VEVs of the field squared and of the energy-momentum tensor which 
diverge on the brane. For Dirichlet boundary condition the current density and its normal derivative 
vanish for points on the brane. The asymptotic behavior of the current density along It h compact 
dimension for large values of the corresponding length Li crucially depends whether the parameter 
defined by (|3.24j) . is zero or not. For / 0 the current density decays exponentially like 

e Llk (i- 1 ) . In the case = 0 the decay is power law, as l/L^ +2u+2 , for both massless and massive 

fields. This behavior for massive fields is essentially different from that for Minkowski bulk where the 
current is suppressed exponentially, by the factor e~ mLl . The expression for the current density in the 
generalized Randall-Sundrum 1-brane model with compact dimensions is obtained from the formulas 
in section [3] with an additional factor 1/2 and with the Robin coefficient (I3.32P for untwisted fields 
and with /3 = 0 for twisted fields. 

The current density in the L-region is given by the expression (13.331) . For a conformally coupled 
massless field, this expression is reduced to the one in Minkowski spacetime with two parallel plates, 
multiplied by the conformal factor ( z/a ) D+1 , with Dirichlet boundary condition on the left plate and 
Robin condition on the right one. On the AdS boundary the brane-induced contribution vanishes as 
z D+ 2 u+ 2 ' p or a observation point, when the location of the brane tends to the AdS horizon, the 

—2znk^ 

brane-induced effects are suppressed by the factor e (i). A similar behavior is exhibited by the 
boundary-free part. From an alternative representation (13.371) it follows that the current density on the 
brane is finite and vanishes for Dirichlet boundary condition. For large values of the length L/, unlike 
the R-region, the decay of the current density in the L-region is exponential for both cases = 0 

and / 0. This feature is related to the discreteness of the spectrum for A. For the value (I2.47|) of 

the Robin coefficient the contribution (13.391) from the special mode has to be added to the right-hand 
side of (|3.37j) . For large values of Li, this contribution dominates in the VEV of the total current 
density. In particular, in the model with a single compact dimension its decay, as a function of Li, is 
power law, as 1/L^ 1 . For the brane location near the AdS boundary, the asymptotic of the current 
density is given by (|3.42l) and it is suppressed by the factor e~' yiLl ^ z °. In the presence of the special 
mode, for a given value of z/zq (fixed distance from the brane), the corresponding contribution to the 
current density behaves as (zq/Li) d ~ 2 and, hence, it dominates in the total VEV. 

The numerical results have been presented for the D = 4 model with a single compact dimension 
and for a minimally coupled massless scalar field. These results show that, depending on the value of 
the Robin coefficient, the presence of the brane can either increase or decrease the current density. In 
particular, in the example considered, the modulus of the current density takes its minimal value for 
Dirichlet boundary condition. 
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A Other representations of the two-point function 


Here we provide representations for the Hadamard function convenient in the investigation of near 
brane asymptoics for the VEV of the current density. First we consider the R-region. In the represen¬ 
tation (| 2 . 22 l) we separate the series over n/ and for the summation use the Abel-Plana type formula 
[2h(34j 

o_ 00 roo 

7 - ^2 g(k)f(\k\)= du[g(u) + g(-u)]f(u) 

1 ni=— 00 

+i [ du[f(iu)-f(-iu)] J2 e uLM_ V ( AJ ) 

8 =± 1 


with k[ defined in (12.151) (formula (IA.1D is reduced to the standard Abel-Plana formula in the special 
case g(x) = 1, oq = 0). After the application (jA.ll) with g(ki) = e lklAx ‘, the Hadamard function is 
presented in the form 

G(x,x') = G RP +2 xT q~i (x, x') + Gi(x, x'). (A.2) 

Here the part Gr p + 2 xTq -i(x, x') comes from the first term in the right-hand side of (lA.lj) and is 
the Hadamard function for the geometry with a single brane in (D + l)-dimensional AdS spacetime 
with spatial topology R v+2 x T q ~ l for which the All dimension is decompactihed. The lengths of the 
remaining compact dimensions are the same: (L p+ \, ..., L;_i, L/ + 1 ,..., Lp- 1 ). The second term on 
the right of (IA.2I) is induced by the compactihcation of the Ah dimension and is given by the expression 
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(A.3) 


where n g _i = (n p +i,..., rty_i, n; + i,..., no-i), k^ 2 = k 2 — k 2 , and the summation over r in the 
exponent goes over r = 1,... ,D — 1, r ^ In deriving this result, we have used the relation 

p—suAx 1 00 

E _ ! = 2 E e ~ nULl ^sh(uAx l + inai). (A.4) 

s==bl n=l 


Note that the part G RP +i xTg -i(x, x') does not contribute to the VEV of the current density along Ah 
dimension. The expression (IA.3I) gives the contribution to the Hadamard function from the modes 
with continuous A. In the presence of a bound state, a similar representation can be found for the 
corresponding contribution. 

Another useful representation for the Hadamard function in the R-region is obtained by using the 
identity (12.231) for the integrand in (1A.3I) . The part in Gi(x,x') coming from the first term in the 
right-hand side of (I2.23P gives the corresponding function in the boundary-free AdS spacetime (the 
brane is absent), denoted here by Gqi(x,x'). In the part induced by the brane, coming from the last 
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term in (|2.23|) . we rotate the integration contour by the angle vr/2 for the term with j = 1 and by the 
angle —7r/2 for the term with j = 2. In this way we get 


Gi(x,x') = G 0 i{x,x')- 
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(A.5) 


Now, we consider the L-region. The corresponding expression for the Hadamard function is given 
by (12.431) . In a way similar to that for the R-region, applying the formula (IA.1I) . for the part induced 
by the compactification of the Zth dimension we find the following representation 


Gi(x,x') = 
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(A.6) 


In this expression, for the summation over n we can use the formula (12.441) . The part coming from 
the first term in the right-hand side of ([2.44)) corresponds to the boundary free contribution and, as 
a result, we obtain the decomposition similar to (IA.5I) where now the brane-induced contribution is 
obtained from that in (IA.5D by the replacements I u {x) K v {x). 
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